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Abstract. A novel numerical method for solving inverse scattering problem with fixed-energy 
data is proposed. The method contains a new important concept: the stability index of the 
inversion problem. This is a number, computed from the data, which shows how stable the 
inversion is. If this index is small, then the inversion provides a set of potentials which differ so 
little, that practically one can represent this set by one potential. If this index is larger than 
some threshold, then practically one concludes that with the given data the inversion is unstable 
and the potential cannot be identified uniquely from the data. Inversion of the fixed-energy 
phase shifts for several model potentials is considered. The results show practical efficiency of 
the proposed method. The method is of general nature and is applicable to a very wide variety 
of the inverse problems. 



1. Introduction 

Let q{x), a; G K^, be a real- valued potential with compact support. Let i? > be a number such 
that q{x) = for |a;| > R. We also assume that q e L^{Bii) , i?^ = {x : \x\ < R, x e M^}. Let 5^ 
be the unit sphere, and a G S^. For a given energy fc > the scattering solution ^fj{x, a) is defined 
as the solution of 



(1.1) Ai! + k^iP - qix)^; = , xeW^ 

satisfying the following asymptotic condition at infinity: 



(1.2) 



ika-x 



(1.3) 

It can be shown, that 



lim 

r — *oo 



|a:|— r 



dv 
dr 



ikv 



ds = 0. 



(1.4) "ipi^i ct) — "ipo + A{a' , a, fc) h o ( - I , as r ^ oo , — = a' r -.^ \x\. 

r \r J r 

The function A(a',Q;, fc) is called the scattering amplitude, a and a' are the directions of the 
incident and scattered waves, and fc^ is the energy, see [|l^, |l2| . 

For spherically symmetric scatterers q{x) = q{r) the scattering amplitude satisfies A(a' , a, fc) ~ 
A{a' -a, fc). The converse is established in [Q. Following ||l^, the scattering amplitude for q = q{r) 
can be written as 
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oo / 

(1.5) A{a',a,k) = Y. Mmra{a' WUa) , 

1=0 m.= -l 

where Yim are the spherical harmonics, normahzed in L'^{S'^), and the bar denotes the complex 
conjugate. 

The fixed-energy phase shifts — tt < Si < tt {6i ~ S{1, k), k > is fixed) are related to Ai{k) (see 
e.g., |l9|) by the formula: 

(1.6) Mk) = ^e''' sm{6i) . 

In section 2 we give, following [1], formulas for calculating fixed-energy phase shifts for piecewise- 
constant compactly supported potentials, let us denote this class of potentials by PC. Since an 
arbitrary integrable potential can be approximated with the prescribed accuracy by a PC potential, 
the class PC is sufliciently large for practical purposes. 

In sections 3 and 4 a novel minimization method, the stability index method, is described. Our 
inversion procedure is based on this method. An important novel feature of this method, which 
seems not have been present in other methods, is the concept of the stability index, which is a 
number characterizing the stability of the numerical inversion. 

Several parameter-fitting procedures were proposed for calculating the potentials from the fixed- 
energy phase shifts, (by Fiedeldey, Lipperheide, Hooshyar and Razavy, loannides and Mackintosh, 
Newton, Sabatier, May and Scheid, Ramm, and others). These works are referenced and their 
results are described in [5] and [10]. Recent works [6]-[9] and [18]-[20] present new numerical 
methods for solving this problem. 

In section 5 numerical examples of the inversion of the fixed-energy shifts are given for three 
potentials. Physical motivation for the choice of these potentials is given and directions for future 
research are suggested. 

Section 6 contains brief conclusions. 



2. Phase Shifts for Piecewise-Constant Potentials 

Phase shifts for a spherically symmetric potential can be computed by a variety of methods, e.g. 
by a variable phase method described in [Q. The computation involves solving a nonlinear ODE 
for each phase shift. However, if the potential is compactly supported and piecewise-constant, then 
a much simpler method described in can be used. It is summarized below. Since the set of 
compactly supported and piecewise-constant potentials is dense in the set L"'^(0, oo) potentials, it 
is quite reasonable to look for an approximate solution to the inverse scattering problem in the 
class of piecewise-constant compactly supported potentials. 

Consider a finite set of points = tq < ri < 7'2 <.-.< r^r = i? and a piecewise-constant 
potential 

(2.1) q{r) — qi, on [r^-i, r^) for i = 1, . . . , N, and g = for r > i?. 

Denote 



(2.2) 



Kj :— k qi , 



STABILITY INDEX METHOD 



3 



where i = 1, . . . , A^, and k is some fixed positive number. Consider the following problem for the 
radial Schrodinger equation: 

(2.3) f^+(^e-lil±}iy,=q^,, hm(2; + l)!!r-'-Vz(r-) = l, 

which we rewrite as: 

(2-4) -^^i^^-^h-' 



on the interval r^-i < r < r^. On [ri_i,ri) one has the following general solution of (2.4) 



(2.5) ipi{r) = A^ji{K,r) + B^ni{)iir), 

where 



(2.6) Jiikr) = y ^J,^,^,{kr) , n,{kr) = sj ^N,+,/,{kr) 

and Ji , Ni are the Bessel and Neumann functions. 



If k| = fc^ — (7i < for some z, then the solution (pi{r) in (2.4) can be expressed through some 
powers of r (for Ki = 0) or the modified Bessel and Neumann functions, and our approach is still 
valid with the appropriate changes. 

From the regularity of (pi at zero one gets Bi = 0. Denote Xi — Bi/Ai, then xi = 0. We are 
looking for a continuously differentiable solution ipi. Thus, the following interface conditions hold: 

Aiji{Kin) + BiniiKiTi) = Ai+iji[Ki+iri) + Bi+ini{Ki+in), 

(2.7) 

Tr^iAij'iiKiri) + Bin[{Kin)] = A,+ij'i{K.i+iri) + B,+in[{Ki+in). 

Therefore 



(2.8) 



Ai+i\ 1 /all al2^ f A. 



Bi+i) Ki+i \ali a22/ \B 



where the entries of the matrix a' can be written explicitly (see ^ for details) . 
Thus 

(2.9) ..,i = 4i±4^, xr.= ^ 

The phase shift 6{k,l) is defined by 

(2.10) ^i{r) ^ ffi^ sin(fcr - ^ + S{k, I)) r ^ oo , 
where F{k,l) is the Jost function. For r > R one has: 

(2.11) (pi{r) - AN+iJi{kr) + BN+ini{kr). 

From ( ^.11 ) and the asymptotics ji{kr) ~ sin(A;r — ni{kr) ~ — cos(/cr — Z7r/2), r oo, 

one gets: 

(2.12) tan(5(fc, = -^^^ = -xn+i ■ 

An+1 
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Finally, the phase shifts of the potential q{r) are calculated by the formula: 
(2.13) 5(fc, Z) = — arctanxjv+i. 

Let qo{r) be a spherically symmetric piecewise-constant potential. Let {6{k, 0};^i be the set of 
its phase shifts for a fixed fc > and a sufficiently large N. Let q{r) be another potential, and let 
{S{k, be the set of its phase shifts. 

The best fit to data function <i>(q, k) is defined by 



(2.14, *(,,,) .siiwyMSiiie. 

The phase shifts are known to decay rapidly with I, see [Q. Thus, for sufficiently large N, the 
function $ is practically the same as the one which would use all the shifts in ( ^.14 ). The inverse 
problem of the reconstruction of the potential from its fixed-energy phase shifts is reduced to the 
minimization of the objective function <I> over an appropriate admissible set. 



3. Stability Index Minimization Method 
Let the minimization problem be 



(3.1) min{$((7) : q e Aadm} 

Let go be its global minimizer. Typically, the structure of the objective function $ is quite 
complicated: this function may have many local minima. Moreover, the objective function in a 
neighborhood of minima can be nearly flat resulting in large minimizing sets defined by 



(3.2) S,^{qe Aad,n ■■ <^{q) < $(9o) + e} 

for an e > 0. 

Given an e > 0, let be the diameter of the minimizing set S^, which we call the Stability 
Index D of the minimization problem ( 3T ) . The usage of the letter D for this index is explained 
in formula (3.6) below. 

One would expect to obtain stable identification for minimization problems with small stability 
indices. However, the minimization problems with large stability indices have distinct minimizers 
with practically the same values of the objective function. If no additional information is known, 
one has an uncertainty of the minimizer's choice. The stability index provides a quantitative 
measure of this uncertainty or instability of the minimization. 

The basic idea of the Stability Index minimization method is to iteratively estimate normalized 
stability indices of a minimization problem, and, based on this information, to conclude if the 
method has achieved a stable minimum. 

A particular implementation of the Stability Index method used here employs a Hybrid Stochastic- 
Deterministic (HSD) approach. The stochastic part explores the entire admissible set, while the 
deterministic local minimization finds the best fit in a neighborhood of the chosen in the stochastic 
part of the search initial guesses. The HSD approach has proved to be successful for a variety of 
problems in inverse quantum scattering (see |l^) as well as in other applications (see §). A 
somewhat different implementation of the Stability Index Method is described in 

We seek the potentials q{r) in the class of piecewise-constant, spherically symmetric real-valued 
functions. Let the admissible set be 



(3.3) Aadm C {{ri,r2, ■ ■ ■ ,rM,qi,q2, . . . ,qM) ■ < n < R , qiow < qm < qhigh} , 
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where the bounds qiow and qhigh for the potentials, as well as the bound M on the expected 
number of layers are assumed to be known. 

A configuration (ri, r2, . . . , tm, 91, 92, • ■ ■ , to) corresponds to the potential 

(3.4) q{r) = q„i , for r„_i < r < r„i , 1 <m < M , 

where rQ — and q(r) = for r > vm = R- 

Note, that the admissible configurations must also satisfy 

(3.5) ri < r2 < ra < • ■ • < TA/ . 

First we describe the global (stochastic) part of the algorithm, which can be called the Iterative 
Reduced Random Search (IRRS) method. This description is followed by its iterative version, and 
a Local Minimization Method (LMM) incorporating a Reduction procedure. 

Let a batch H oi L trial points be generated in Aadm using a uniformly distributed random 
variable. In our case Aadm. is a box in . The uniform random variable is called 2M times 
to produce a point (configuration representing a potential) in this box (after the appropriate 
rescaling in each dimension). Finally, the obtained values of are rearranged in the ascending 
order to satisfy ( |3.5|) . 

A certain fixed fraction 7 of the original batch of L points is used to proceed with the local 
searches. Typically, L = 5000 and 7 = 0.01. This reduced sample Hred of 7L points is chosen to 
contain the points with the smallest jL values of $ among the original batch H. The local searches 
(the LMM procedure) are started from every point in this reduced sample Hred ■ This way only 
the points that seem to be in a neighborhood of the global minimum are used for an expensive 
local minimization, and the computational time is not wasted on less promising candidates. 

Let Hmin be the 7L points obtained as the result of the local minimizations (7L = 50 in our 
computations). Let Smin be the subset of Hmin containing points {pi\ with the smallest v^L 
(0 < < 1, we used v — 0.16) values in Hmin- We call Smin the minimizing set. The choice of u 
determines a representative sample of global minimizers. If all these minimizers are close to each 
other, then the objective function <i> is not flat near the global minimum. That is, the method 
identifies the minimum consistently. To define this consistency in quantitative terms, let ||.|| be a 
norm in the admissible set. 

Let 

(3.6) D = diam{Smin)^'aia2i{\\Pi-Pj\\/dav ■ Pi.Pj^Smm}, 

where dav is the average norm of the elements in Hmin- The normalization by dav is introduced 
to provide comparable results for different potentials. Thus D is an estimate for the (normalized) 
Stability Index of the minimization problem. The identification is considered to be stable if the 
Stability Index D < e. Otherwise, another batch of trial points is generated, and the process is 
repeated as follows. 

Iterative Reduced Random Search (IRRS). (at the j— th iteration). 
Fix <j,iy < I, (3 > I, e>0 and jmax- 

1. Generate another batch H^ of L trial points (configurations) in Aadm using a uniform random 
distribution. 

2. Reduce H^ to the reduced sample Hij,^^ of 7L points by selecting the points in H^ with the 
smallest 7L values of <i>. 

3. Apply the Local Minimization Method (LMM) starting it at each of the jL points in iJ^g^, 
and obtain the set H^^^ consisting of the 7L minimizers. 
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4. Combine i?^„i„ with H^^^ obtained at the previous iteration. Let S'^i„ be the set of v^L 
points from U H^^^-^^ with the smallest values of $. (Use i?,^„j„ for j = 1). 

5. Compute the Stability Index (diameter) of >S';^j„ by = max{||pi —pk\\/dav '■ Pi,Pk € 
Smin} , where dav = J2 \\Pk\\/lL , pk e H^in is the average norm of the elements of i?^jj„. 
(Thus, dav is computed only once at j = 1 and this value is used for all subsequent iterations). 

6. Stopping criterion. 

Let p e ■S'^^jjn be the point with the smallest value of $ in S'^jj„ (the global minimizer) . 
If < e, then stop. The global minimum is p. The minimization is stable. 
If > e and $((7) < (3^{p) : q G *5'/„i„7 then stop. The minimization is unstable. The 
Stability Index is the measure of the instability of the minimization. 

Otherwise, return to step 1 and do another iteration, unless the maximum number of 
iterations jmax is exceeded. 
We used (3 = 1.1, e = 0.02 and jmax = 30. The choice of these and other parameters (L = 
5000, 7 = 0.01, 1/ — 0.16 Er = 0.1 (used in LMM)) is dictated by their meaning in the algorithm and 
the comparative performance of the program at their different values. As usual, some adjustment 
of the parameters, stopping criteria, etc., is needed to achieve the optimal performance of the 
algorithm. 

4. Local Minimization Method 

The Hybrid Stochastic Deterministic Method couples the Stochastic part described in the pre- 
vious section with a deterministic Local Minimization Method. Numerical experiments show that 
the objective function $ is relatively well behaved in this problem: while it contains many local 
minima and, at some points, $ is not differentiable, standard minimization methods work well 
here. A Newton-type method for the minimization of $ is described in [Q. We have chosen to use 
a variation of Powell's minimization method which does not require the computation of the deriva- 
tives of the objective function. Such method needs a minimization routine for a one-dimensional 
minimization of <&, which we do using a Bisection or a Golden Rule method. See or for a 
complete description of our method. 

Now we can describe our Basic Local Minimization Method in E^*^, which is a modification 
of Powell's minimization method [|j . It is assumed here that the starting position (configuration) 
Qo E Aadm is suppled by the procedure LMM (see below), and the entry to LMM is provided by 
the global minimization part (IRRS). 

Basic Local Minimization Method. 

1. Choose the set of directions ui , i = 1, 2, . . . , 2M, to be the standard basis in 9?^-^ 

u, = (0,0,...,l,...,0), 

where 1 is in the i-th place. 

2. Save your starting configuration supplied by LMM as Qo • 

3. For each i = 1, . . . , 2M move from Qo along the line defined by Ui and find the point of 
minimum Q\. This defines 2M temporary points of minima. 

4. Re-index the directions Ui, so that (for the new indices) <&(Qi) < ^{Q^) <i ■ • ■ i '^{Q2m) — 

HQo)- 

5. For i = 1, . . . , 2M move from Qi-i along the direction Ui and find the point of minimum Qi. 

6. Set V = Q2M - Qo- 

7. Move from Qo along the direction v and find the minimum. Call it Qo again. It replaces Qo 
from step 2. 

8. Repeat the above steps until a stopping criterion is satisfied. 
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Note, that we use the temporary points of minima Q\ only to rearrange the initial directions Ui 
in a different order. The stopping criterion is the same as the one in ||ll|, Subroutine Powell]. 
Still another refinement of the local phase is necessary to produce a successful minimization. 



The admissible set Aadm, see (5.3)-(3.5), belongs to a 2M dimensional minimization space M^^^. 
The dimension 2M of this space is chosen a priori to be larger than 2N , where N is the number 
of layers in the original potential. We have chosen M = 2 in our numerical experiments. However, 
since the sought potential may have fewer than M layers, we found that conducting searches in 
lower-dimensional subspaces of is essential for the local minimization phase. A variation of 
the following "reduction" procedure has also been found to be necessary in [|| for the search of 
small subsurface objects, and in |^ for the identification of multilayered scatterers. 

If two adjacent layers in a potential have values Vi-i and Vi and the objective function $ is 
not changed much when both layers are assigned the same value Vi (or then these two 

layers can be replaced with just one layer occupying their place. The change in <I> is controlled 
by the parameter e,-. We used = 0.1. This value, found from numerical experiments, seems to 
provide the most consistent identification. The minimization problem becomes constrained to a 
lower dimensional subspace of and the local minimization is done in this subspace. 

Reduction Procedure. Let be a positive number. 

1. Save your starting configuration Qo ~ {ri,r2, ■ ■ ■ ,rM ,vi,V2, ■ ■ ■ , vm) G Aadm and the value 
$(Qo)- Let the (M + l)-st layer be Lm+i = {"Tm < < R} and um+i — 0. 

2. Let 2 < z < M + 1. Replace Vi-i in the layer L^-i by Vi. This defines a new configuration 
Q^, where the layers Li-i and Li are replaced with one new layer. Here d stands for the 
downward adjustment. Compute and the difference cf = |$(Qo) ^ ^{Qi)\- Repeat 
for each layer in the original configuration Qo- 

3. Let 1 < i < M . Replace Vi^i in the layer Li+i by Vi. This defines a new configuration Q", 
where the layers Li and Li+i are replaced with one new layer. Here u stands for the upward 
adjustment. Compute and the difference c" = |<i>((5o) — ^{Qi)\- Repeat for each layer 
in the original configuration Qo- 

4. Find the smallest among the numbers cf and c". If this number is less than (^^{Qq), then 
implement the adjustment that produced this number. The resulting new configuration has 
one less layer than the original configuration Qo- 

5. Repeat the above steps until no further reduction in the number of layers is occurring. 
Note, that an application of the Reduction Procedure may or may not result in the actual 

reduction of the number of layers. 

Finally, the entire Local Minimization Method (LMM) consists of the following: 

Local Minimization Method (LMM). 

1. Let your starting configuration supplied by IRRS be Qq — (ri, r2, . . . , tm, wi, W2, • ■ ■ , vm) G 

Aadm • 

2. Apply the Reduction Procedure to Qo, and obtain a reduced configuration Qq containing 

layers. 

3. Apply the Basic Minimization Method in AadmCl^'^^^ with the starting configuration Qq, 
and obtain a configuration Qi. 

4. Apply the Reduction Procedure to Qi, and obtain a final reduced configuration Q\. 

As we have already mentioned, LMM is used as the local phase of the global minimization. 

5. Numerical Results 

We studied the performance of the algorithm for 3 different potentials qi{r), i — 1,2,3 chosen 
from the physical considerations. 
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The potential 93 (r) = — 10 for < r < 8.0 and qs = for r > 8.0 and a wave number k — 1 
constitute a typical example for elastic scattering of neutral particles in nuclear and atomic physics. 
In nuclear physics one measures the length in units of fm = 10~^^m, the quantity in units of 
1/fm^, and the wave number in units of 1/fm. The physical potential and incident energy are 
given by V{r) = |^'?3('") and E ~ respectively, here h :— h — 6.62510"^'' erg-s is the 

Planck constant, he = 197.32 MeV-fm, c = 3 ■ 10^ m/sec is the velocity of light, and fi is the mass 
of a neutron. By choosing the mass ^ to be equal to the mass of a neutron /i = 939.6 MeV/c^, the 
potential and energy have the values of V{r) = -207.2 MeV for < ?■ < 8.0 fm and E{k =l/fm ) 
= 20.72 MeV. In atomic physics one uses atomic units with the Bohr radius Oq = 0.529 • lO^^^m 
as the unit of length. Here, r, k and q^ are measured in units of oq, 1/ao and l/og, respectively. 
By assuming a scattering of an electron with mass mg = 0.511 MeV/c^, we obtain the potential 
and energy as follows: V{r) = -136 eV for < r < 8ao = 4.23 • 10-^°m and E{k = 1/ao) = 13.6 
eV. These numbers give motivation for the choice of examples applicable in nuclear and atomic 
physics. 

The method used in this paper deals with finite-range (compactly supported) potentials. One 
can use this method for potentials with the Coulomb tail or other potentials of interest in physics, 
which are not of finite range. This is done by using the phase shifts transformation method which 
allows one to transform the phase shifts corresponding to a potential, not of finite range, whose 
behavior is known for r > a, where a is some radius, into the phase shifts corresponding to a 
potential of finite range a (see [2], p. 156). 

In practice differential cross section is measured at various angles, and from it the fixed-energy 
phase shifts are calculated by a parameter-fitting procedure. Therefore, we plan in the future work 
to generalize the stability index method to the case when the original data are the values of the 
differential cross section, rather than the phase shifts. 

By the physical reasons discussed above, we choose the following three potentials: 



91 W 



-2/3 < r < 8.0 
0.0 r > 8.0 



-4.0 < r < 8.0 
0.0 r>8.0 



-10.0 < r < 8.0 
0.0 r>8.0 



In each case the following values of the parameters have been used. The radius R of the support 
of each was chosen to be i? = 10.0. The admissible set Aadm ( |3.3| ) was defined with AI — 2. The 
Reduced Random Search parameters: L — 5000 , 7 = 0.01 , v — 0.16 , e = 0.02 , /? — 1.10 ,jmax — 
30. The value = 0.1 was used in the Reduction Procedure during the local minimization phase. 
The initial configurations were generated using a random number generator with seeds determined 
by the system time. A typical run time was about 10 minutes on a 333 MHz PC, depending on the 
number of iterations in IRRS. The number TV of the shifts used in ( 2.14| ) for the formation of the 



objective function ^{q) was 31 for all the wave numbers. As it can be seen from Table 1 the shifts 
for the potential 93 decay rapidly for fc = 1, but they remain large for fc = 4. The upper and lower 
bounds for the potentials qiow = —20.0 and qhigh = 0.0 used in the definition of the admissible set 
■Aadm were chosen to reflect a priori information about the potentials. 
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Table 1. Phase shifts 6{k,l) of qsir). 



I 




= 1.0 


k 


= 4.0 





-0. 


.66496 


-0, 


.62217 


1 


-0, 


,31009 


-0, 


.64598 


2 


-0. 


,72324 


-0, 


.65824 


3 


-0. 


88586 


-0, 


.64604 


4 


-0, 


,74713 


-0, 


.74239 


5 


1, 


.15839 


-0, 


65260 


6 


1, 


.54292 


-0, 


86826 


7 


-0, 


,56945 


-0, 


.69851 


8 


-0, 


,38745 


_1, 


00663 


9 


0, 


.16888 


-0, 


.83757 


10 


-0, 


02690 


_1, 


.08641 


11 


-0, 


,01261 


_1, 


.09074 


12 


0, 


.00017 


-1, 


.08645 


13 


-0, 


,00010 


_1, 


39603 


14 


-0, 


,00004 


-1, 


.24536 


15 


0, 


oofino 


_1 


55399 


16 


0, 


ooooo 


1, 


,49036 


17 


0, 


ooooo 


1, 


,56437 


18 


0, 


OOOOO 


1, 


,11836 


19 


0, 


,00000 


1, 


,12265 


20 


0, 


,00000 


1, 


,11829 


21 


0, 


,00000 


0, 


,60125 


22 


0, 


,00000 


0, 


,58327 


23 


0, 


,00000 


0, 


,59973 


24 


0, 


,00000 


0, 


,00875 


25 


0, 


.00000 


-0, 


.18826 


26 


0, 


,00000 


-0, 


.11221 


27 


0, 


,00000 


-0, 


.47503 


28 


0, 


.00000 


-1, 


.22725 


29 


0, 


.00000 


-1, 


.29222 


30 


0, 


.00000 


-1, 


.25626 



The identification was attempted with 3 different noise levels h. The levels are h — 0.00 (no 
noise), h = 0.01 and h = 0.1. More precisely, the noisy phase shifts 6h{k,l) were obtained from 
the exact phase shifts 6{k, I) by the formula 

dh{k,l) = d{k,l)il + {0.5 - z) ■ h) , 

where z is the uniformly distributed on [0, 1] random variable. 

The distance d{pi{r),p2{r)) for potentials in step 5 of the IRRS algorithm was computed as 

d{pi{r),p2{r)) = \\pi{r) - P2{r)\\ 

where the norm is the L2-norm in . 

The results of the identification algorithm (the Stability Indices) for different iterations of the 
IRRS algorithm are shown in Tables 2-4. 
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Table 2. Stability Indices for qi{r) identification at different noise levels h. 



k 


Iteration 


h = 0.00 


h = 0.01 


h = 0.10 


1.00 


1 


1.256985 


0.592597 


1.953778 




2 


0.538440 


0.133685 


0.799142 




3 


0.538253 


0.007360 


0.596742 




4 


0.014616 




0.123247 




5 






0.015899 


2.00 


1 


0.000000 


0.020204 


0.009607 


2.r>{) 


1 


o.oooooo 


{).{)lior,:i 


0.046275 


3.00 


1 


0.000000 


0.000501 


0.096444 


4.00 


1 


0.000000 


0.022935 


0.027214 



For example, Table 4 shows that for k = 2.5, h = 0.00 the Stability Index has reached the 
value 0.013621 after 2 iteration. According to the Stopping criterion for IRRS, the program has 
been terminated with the conclusion that the identification was stable. In this case the potential 
identified by the program was 



p{r) 



-10.000024 
0.0 



< r < 7.999994 
r > 7.999994 



which is very close to the original potential 



93 (r) = 



f-10.0 0<r<8.0 
0.0 r>8.0 



On the other hand, when the phase shifts of qs{r) were corrupted by a 10% noise (fc = 2.5, h = 
0.10), the program was terminated (according to the Stopping criterion) after 4 iterations with the 
Stability Index at 0.079241. Since the Stability Index is greater than the a priori chosen threshold 
of e = 0.02 the conclusion is that the identification is unstable. A closer look into this situation 
reveals that the values of the objective function <&(pi), Pi G Smin (there are 8 elements in Smin) 
are between 0.0992806 and 0.100320. Since we chose /3 = 1.1 the values are within the required 
10% of each other. The actual potentials for which the normalized distance is equal to the Stability 
Index 0.079241 are 



-9.997164 
Pi(r) = I -7.487082 
0.0 



< r < 7.932678 
7.932678 < r < 8.025500 
r > 8.025500 



and 



P2{r) 



-9.999565 < r < 7.987208 
-1.236253 7.987208 < r < 8.102628 
0.0 r > 8.102628 



with $(j3i) = 0.0992806 and $(^2) = 0.0997561. One may conclude from this example that the 
threshold e = 0.02 is too tight and can be relaxed, if the above uncertainty is acceptable. 
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Table 3. Stability Indices for q2{r) identification at different noise levels h. 
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k 


Iteration 


h = 0.00 


h = 0.01 


h = 0.10 


1.00 


1 


774376 


598471 


108902 




2 


77371 8 


1 027345 


023206 




3 


0.026492 


0.025593 


0.023206 




4 


0.020522 


0.029533 


0.024081 




5 


0.020524 


0.029533 


0.024081 




6 


0.000745 


0.029533 






7 

8 




029533 

029533 






9 




0.029533 






10 




0.029533 






11 




0.029619 






12 




0.025816 






13 




0.025816 






14 




0.008901 




2.00 


1 


0.863796 


0.799356 


0.981239 




2 


0.861842 


0.799356 


0.029445 




3 


0.008653 


0.000993 


0.029445 




4 






0.029445 




5 






0.026513 




6 






0.026513 




7 






0.024881 


2.50 


1 


1 848910 


1.632298 


894087 




2 


1.197131 


1.632298 


507953 




3 


580361 


1 183455 


025454 




4 


030516 


528979 






5 


016195 


032661 




3.00 


1 


1 844702 


1 849016 


1 708201 




2 


1 649700 


1 782775 


1 512821 




3 


1 456026 


1 782775 


1 412345 




4 


1 410253 


1 457020 


1 156964 




5 


624358 


961263 


1 156964 




6 


692080 


961263 


902681 




7 


692080 


961263 


902681 






o 




n 9Q1 fil 1 


u.yuzft / ft 




9 


0.345804 


0.286390 


0.159221 




10 


0.345804 


0.260693 


0.154829 




11 


0.043845 


0.260693 


0.154829 




12 


0.043845 


0.260693 


0.135537 




13 


0.043845 


0.260693 


0.135537 




14 


0.043845 


0.260693 


0.135537 




15 


0.042080 


0.157024 


0.107548 




16 


0.042080 


0.157024 






17 


0.042080 


0.157024 






18 


0.000429 


0.157024 






19 




0.022988 





4.00 



1 0.000000 0.000674 0.050705 
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Table 4. Stability Indices for q^ir) identification at different noise levels h. 



k 


Iteration 


h = 0.00 


h = 0.01 


h = 0.10 


1.00 


1 


564168 


594314 


764340 




2 


0.024441 


028558 


081888 




3 


0.024441 


014468 


050755 




4 


0.024684 








5 


0.024684 








6 


0.005800 






2.00 


1 


684053 


1.450148 


485783 




2 


0.423283 


792431 


078716 

yj »\j i ^ i -1- yj 




3 


006291 


457650 


078716 




4 




0.023157 


0.078716 




5 






0.078716 




6 






0.078716 




7 






0.078716 




8 






0.078716 




9 






0.078716 




10 






0.078716 
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0.078716 


2.50 


1 


0.126528 


0.993192 


0.996519 




2 


0.013621 


0.105537 


0.855049 




3 




0.033694 


0.849123 




4 




0.026811 


0.079241 


3.00 


1 


0.962483 


1.541714 


0.731315 




2 


0.222880 


0.164744 


0.731315 




3 


0.158809 


0.021775 


0.072009 




4 


0.021366 








5 


0.021366 








6 


0.001416 






4.00 


1 


1.714951 


1.413549 


0.788434 




2 


0.033024 


0.075503 


0.024482 




3 


0.018250 


0.029385 






4 




0.029421 






5 




0.029421 






6 




0.015946 





6. Conclusions 

A novel numerical method for solving inverse scattering problem with fixed-energy data is pro- 
posed. The method contains a new important concept: the stability index of the inversion problem. 
This index is a number, computed from the data, which shows how stable the inversion is. If this 
index is small, then the inversion provides a set of potentials which differ so little, that practically 
one can represent this set by one potential. If this index is larger than some threshold, then one 
concludes that practically, with the given data, the inversion is unstable and the potential cannot 
be identified uniquely from these data. Inversion of the fixed-energy phase shifts for several model 
potentials is considered. The results show practical efficiency of the proposed method. The method 
is of general nature and is applicable to a very wide variety of inverse problems. 
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